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ABSTRACT: 

The key element in Mechatronics design is the concurrent synergetic integration, modeling, 
simulation, analysis and optimization of multidisciplinary knowledge through the design process from 
the very start of the design process. Mechatronics engineer is expected to design engineering systems 
with synergy and integration toward constrains like higher performance, speed, precision, efficiency, 
lower costs and functionality. This paper proposes the conception and development of ball and beam 
system based on Mechatronics design approach. A complete overall system and subsystems selection, 
modeling, simulation, analysis, and integration are presented. The proposed Mechatronics design and 
models were created and verified using MATLAB /Simulink software and are intended for research 
purposes, as well as, application in educational process. 
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I. INTRODUCTION 

Mechatronics is defined as synergistic integration of mechanical engineering, electric engineering, 
electronic systems, information technology, intelligent control system, computer hardware and software to 
manage complexity, uncertainty, and communication in design and manufacture of products and processes, a 
working definition can be; the synergistic integration of sensors, actuators, signal conditioning, power 
electronics, decision, control algorithms, computer hardware and software to manage complexity, uncertainty, 
and communication in engineered systems. The key element in Mechatronics design is the concurrent synergetic 
integration (instead of sequential), analysis and optimization of these areas through the design process resulting 
in products with more synergy [1] and a balance between modeling/analysis and experimentation /hardware 
implementation. Based on this, Mechatronics engineer is expected to design engineering systems with synergy 
and integration toward constrains like higher performance, speed, precision, efficiency, lower costs and 
functionality. Reference to [1-2], Mechatronics system design process can be divided into systematic, simple 
and clear design steps including; Problem statement; Conceptual design and functional specifications; Parallel 
(concurrent) design and integration of system and all subsystems as whole and concurrently including; selection, 
design and synergetic integration of mechanical, electronics, software, control unit, control algorithm and 
interface sub-systems; Modeling and simulation; Prototyping, Testing and Optimization; and finally 
Manufacturing, and commercialization, these steps will be followed and explained in this paper, were this paper 
proposes the conception and development of ball and beam system based on Mechatronics design approach. A 
complete overall system's and subsystem's components' selection, design, integration, as well as, modeling, 
simulation and analysis are to be presented. The significance of the ball and beam system is due to it's 
complexity, and at the same time it's simplicity, where it is a double integrator system that is open-loop unstable 
and simply, it is a position control system, presents a challenging design and control problem, it is often used by 
many control theorists and engineers as a bench mark problem for testing and analyzing the results of many 
different theoretical control schemes and new controller and observer methods. 

II. PRE-STUDY PROCESS-PROBLEM STATEMENT 

The ball and beam system is motion control system consisting of three main parts; control system, a 
servo or electric motor , free rolling back and forth, Ball on top of long V-grooved Beam whose inclination can 
be adjusted by manipulating the angular position of the servo. Ball and beam system is often used as a bench 
mark problem for many different control schemes and presents a challenging design and control problem, due to 
fact that it is double integrator (two poles at the origin) system that is open-loop unstable, in such system even at 
optimal beam position conditions the ball will swing on Beam, up to roll off the end of the beam. Target user; 
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Researchers and Educators. User 's requirements; it is required to develop a ball and beam system; linear, one- 
dimensional single-variable experimental unit, that can be used for testing, analysis, education, and 
demonstration the effect of various control strategies, mainly SISO systems, precise, cost-efficient, with simple, 
easy to understand and use interface allowing the users to manipulate a machine, gathering and outputting data, 
space saving, and with attractive design. System parameters, requirements and analysis, are listed in Table 1. 

III. CONCEPTUAL DESIGN 

The purpose of design is to propose the conception and development of ball and beam system based on 
Mechatronics design approach, it is required to propose a design of a ball and beam system that fits user's 
(Researchers, Educators , students) needs and requirements, with overall function ( intended) to help users to 
test , visualize and analyze the selected and designed control strategy to continually measure and stabilize the 
rolling ball position a long a given beam by adjusting the rotating angle of a beam by tilting it frontwards, 
backwards, leftwards, or rightwards, all this accomplished by manipulating the angular position of a selected 
servo. System's functional structure block diagram is shown in Figure 1. Morphological table, analysis and 
evaluating the best solution are shown in Table 2. Block diagram representations, preliminary structures 
including mechanical, electronic and system dynamics are shown in Figure 2 and Figure 3 
Table 1 requirements analysis 



Requirements 


Value/ 
Unit 


Requirements type 


Fixed 
Req. 


Soft 
Req. 


Quantitative 
Req. 


Qualitative 
Req. 


Control system 


Compact, cost effective. 




Fixed 






Qual. 


Control 

;il <JOr[| hill 

UIIlUI 1(11111 


Simple, quick, precise, 
efficient, easy to program 






Soft 






Beam 


V-grooved, 


40 cm 


Fixed 




Quan. 




Mass , m 


1.5 kg 


Beam length, m 


40 cm 


Ball 


Radius 


0.015 m 


Fixed 








Mass, M 


1111™ 


Actuators 


One actuator, for adjusting 
the beam rotating angle, 
available, inexpensive, 
simple to interface, and 
control. 


12-24 v 




Soft 




Qual. 


Sensors 


Position sensor, simple to 
interface, available, 
inexpensive 






Soft 




Qual. 


Transmission 


Gears, belt, chain 
Output gear radius 


0.03 m 




Soft 






Overall gear ratio 


n 


Overall System 
dimensions 


Height: 


40 Cm 


Fixed 






Qual. 


Width 


30 Cm 


Length 


50 Cm 


User Interface 


Simple and easy to use and 
understand 










Qual. 


Machine 

aesthetics 

design 


Compact, attractive fits 
user requirements. 






Soft 


Quan 
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Tunc don 


Solution 1 


Solution 2 


Solution 3 


Soludon4 


Control system 


^ficro c ontrollerjl 


Microcomputer 


Analogue 




Control algorithm: 

Position ccniro! 


PID 


rPD : deadbeat7\ 
S^pre filt er^ii' 


Fuzzv 




Ball : steel 


[Rigid ^-^-^^ 


hofliiB 






Beam: shnnumim 


Two parallel thin( 
bats / 


V- grooved J 


U-grooved 




_ »k_ Llld III E \ 


FK -1 mritriT 

Lf '■- I£1U LUI 1 








Position Sensor 


Potentiometer 


Encoder 


Ultia-,ulJ J 


^Linear t saistive vvii e 


Transmission 


Belt 


Chains 


PI — fry p ihla 


t£ears and shaftj 


Interface ( 


Touch screen 


rue vp act 


Switches 




Support 


Slides 


Linkage -j oint 


\^e els 




Poiver 1 


JHectric _^Jnw 








Operator 


Standing 1 


Seating at front J 


Seating at rear 


Remote control 



Table 2 Morphological table, analysis and evaluating the best solution selection 



Switching system on-off 



User interface 



Control system 
Actuator, electric 
energy 

Gears , lever arm 



Position sensor 



Position sensor 



Desired ball position 



Actuator motion 
response 



Motions transmission 



Beam angular position 



Ball position & 
measuring 



Figure 1 Functional Structure block diagram 



IV. PARALLEL (CONCURRENT) SELECTION, DESIGN AND INTEGRATION OF 

SUB-SYSTEMS AND OVERALL SYSTEM. 

Mechatronics design is the synergistic selection, evaluation and integration of the system and all its 
components as a whole and concurrently including sensors, actuators, signal conditioning, power electronics, 
decision, control algorithm, computer hardware and software. All the design disciplines work in parallel and 
collaboratively throughout the design and development process to produce an overall optimal design. The ball 
and beam system is to be divided into realizable hardware and software subsystems including; Mechanical 
subsystem; sensor subsystem ; actuator subsystem ; power supplies, controller and Control algorithm, drive, and 
conditioning circuits), the optimal selection, modeling, simulation, integration, optimization and the exchange of 
information between different modules are to be designed concurrently. 



4.1 The optimal mechanical design: 

The mechanical design is the skeleton of Mechatronics systems, there are several arrangements of 
constructing ball and beam system, three system's arrangements are shown in Figure 2. In Figure 2 (a), one end 
of the beam is coupled to electric motor through lever arm and gears, the other end is fixed, and supported with 
vertical beam. In Figure 2(c) one end of the beam is coupled to electric motor through lever arm and gears, the 
other end is free to move up and down, the beam is supported with vertical bar at its center. In Figure 2(d) the 
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beam at it center, is mounted on the output shaft of an electric motor. System CAD model is shown in Figure 
2(e). The whole system layout is shown in Figure 2(a), and the mechanical system parts, materials and optimal 
dimensions are listed in Table 3. Preliminary block diagram and layout representations of proposed system and 
main components are shown in Figure 2(f)), 



Lever arm 




Gear radius 



9, output 
jsa: angle 




Output gear angle 



Motor with gears 



Output s 
DC motor & gear 



Figure 2(a) One end of the beam is coupled to electric motor through lever arm and gears, the other end is fixed. 




Figure 2 (b) Electric 
motor with gears 




Figure 2 (c) support at beam's 
center 




Figure 2 (d) The beam at it center, is 
mounted on the output shaft of an electric 
motor. 




Figure 2(e) 3-D CAD model 
D 



User 
interface 



Ball ,VV Actuator UAJ ^ 
Position,! hJ controller controller -/-vy 




a=d/L 




Figure 2 (f) system block diagram 
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Table 3 mechanical system parts, type of joints materials and dimensions 





Value: 


Unit 


System dimensions 


Height 


40 Cm 


Width 


30 Cm 


Length 


50 Cm 


Beam 


Shape & material 


V-grooved , steel 


Beam mass, m 


1.5 kg 


Beam length, L 


40 cm 


Ball 


Shape & material 


Steel 


Radius 


0.015 m 


Mass, M 


0.111 kg 


Transmission 


Outer gear radius 


0.03 m 


Overall gear ratio 


n 


Joints 


Re volute joint 




Support 


Steel bar 




Type of joints 


Revolute joints 


3J 


Base, Frame 







4.2.Controller and control algorithm selection 

Due to fact that Ball and beam system is double integrator (two poles at the origin) system, the 
behavior of moving Ball on a beam is behavior of open loop unstable system, in order to stabilize the double 
integrator and open-loop unstable ball and beam system, a special control system is to be designed. The control 
of unstable systems is critically important to many of most difficult control problems and must be studied in 
laboratory, the problem is that real unstable systems are usually dangerous and can not brought in laboratory, 
such unstable systems include; the control model in the rocket or aircraft toppling control system during vertical 
take-off, it is the control where a feedback system is used to prevent rocket or aircraft to topple out of balance 
during launch by forces and moments that could perturb the vertical motion, the angle of thruster jets or 
diverters must be continually controlled to prevent the rocket tumbling or aircraft tipping [2-4]. The ball and 
beam system was developed to resolve this paradox, it is simple, safe mechanism and yet it has the important 
dynamic features of unstable system [2],. 

Control algorithm selection is the most critical decision in the Mechatronics design process, there are number 
of alternative control algorithm strategies that can be applied to the ball and beam system, including but not 
limited to ; PD control, lead compensator, Fuzzy control, Robust control, Linear quadratic Gaussian, Linear 
quadratic regulator. The control system task, in ball and beam system, is to continually measure and stabilize 
the rolling ball position a long the beam by adjusting the rotating angle of a beam ( tilting it frontwards or 
backwards) by manipulating the angular position of the servo, in such a way that the actual position of the ball 
reaches desired position, this is difficult control task, because even at optimal beam position conditions the ball 
will swing on Beam with acceleration proportional to tilt angle, this behavior is behavior of open loop unstable 
system, where the system output (ball position on beam) increases without limits for fixed input (beam angle) 
and is a challenging design and control problem. In [5] W. Yu and F. Ortiz studied complete nonlinear model of 
beam and ball system and stability analysis of PD control. In [6] a designed and implemented both conventional 
pole-placement and neural network methods for beam and ball system are introduced. In [7] the robust nonlinear 
servomechanism theory was applied to design a tracking controller for beam and ball system. In [8] switching 
control for beam and ball system is studied. In [9] fuzzy logic controller and real time implementation of beam 
and ball balancing system is studied. In this paper will be introduce PD controller with deadbeat response, PID, 
and Lead integral compensator. 

Controller selection :A vast variety of controllers (physical controller) are available in the market, embedded 
Microcontroller is optimal selection, since it is inexpensive single chip computer, easy to embed into larger 
electronic circuit designs, also, because of their versatility, Microcontrollers add a lot of power, control, and 
options at little cost; capable of storing and running programs, programmed to perform a wide range of control 
tasks. Optimal microcontroller is PICmicro Microcontroller, supplied with 5VDC. 
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4.3 Actuator (electric machine) selection and integration 

Actuator converts an information signal from the microcontroller control unit, into energy acting on the 
basic system. Mechatronics systems often use electric motors to drive their work loads, the electric actuator 
most suitable to be used for ball and beam is DC motor, due to simple principle of working, quick instantaneous 
and accurate torque generation, available, inexpensive, reversible, and ease designing and implementing 
controller to achieve optimal instantaneous, precise motion control performance. PMDC motor turns electrical 
energy into mechanical energy and produces the torque required to move (rotate) the beam to the desired 
position, L , 

Actuator placement and integration: to physically integrate Mechatronics system components, and to give the 
overall system safe, compact and attractive design , the selected DC motor and gears, are embedded within 
mechanical design and to be located on adjacent ( right side, a whole system base) cuboide housing, with slots 
for digital user interface and input/output ports(see Figure 3) . 



Revolute joint 



r 




Resistive wire sensor DC motor, gears 



Lever arm 




Resistive wire sensor 



V-grooved beam 



User interface EEI 
II 



Electronic circuitry 



DD 



Figure 3(a) 



Figure 3(b) 



Revolute joint 



Resistive wire sensor 



Lever arm 




Figure 3(c) 

Figure 3(a)(b)(c) Three arrangements of overall system design, components and integration 
4.4 Sensor selection 

Sensor converts a state variable of the basic system, into an information signal to the control unit. To 
calculate the error between the actual and desired positions, we need to convert the actual position, into voltage, 
V, and then compare this voltage with the input voltage; the difference between both is the error signal in volts. 
Ball and beam system, provides two outputs measurements and correspondingly two sensors, one to measure 
ball translational position, x a long the beam and other to measure motor shaft angular position 0, 
correspondingly, the optimal control system design requires two control loops designs, inner for motor angle 
control, to make motor follow the input angle and other outer loop to control the ball position on the beam, to 
make ball position follow the reference input position. The two angular quantities the beam angle, a and DC 
motor output position angle, 9 are not the same, but are correlated in simple relationship, relating these two 
quantities result in that only one sensor can be used to measure the position. The position of the ball on the beam 
can be measured using a special sensor, e.g. linear Potentiometer, infrared, and a linear resistive wire, the 
optimal selection is linear resistive wire, to used as a position sensor and the conductive ball moving along it 
acting as a voltage divider ( see Figure 3). 
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4.5 Output signal, Conditioning and Interfacing selection, design and integration 

A most suitable and simple to implement drives for PMDC motor are H-bridge and IC circuit 
e.g.L93D.The H-bridge circuit is supplied with 12VDC and the four bits outputs of microconttoller to drive the 
desire conditions of electric Motor. By H-bridge four NPN Power transistors are used as switch to control the 
direction of current flows to the Motor, and correspondingly conttol the angular position of beam, selected 
microcontroller type is supported with ADC pins, to convert the analog input of sensor readings to digital value, 
fed to microcontroller. A common carrier (see Figure 3-4) for electronics, data outputting and control circuitry 
to be located inside the cuboide housing, that will include slots for user interface and input/output ports 



«- User interface 



Sensor : ball position 



d/L i 



5V 

o 



Figure 4 Microcontroller based control Circuitry and interfacing diagram 

V. Modeling, simulation, analysis and evaluation 

The key essential characteristics of a Mechatronics engineer are a balance between modeling/analysis 
skills and experimentation/hardware implementation skills. Modeling, simulation, analysis and evaluation 
processes in Mechatronics design consists of two levels, sub-systems models and whole system model with 
various sub-system models interacting similar to real situation, all engineering subsystems should be included in 
whole system model. Based on the specification of requirements and design, the subsystems models and the 
whole system model, are to be tested and analyzed, specifications to test and check whether the given design 
specifications are satisfied, If the specification are not satisfied,, modifications can be made, if the specifications 
are satisfied the model can be Optimized. The derivation of complete and accurate ball and beam system's 
dynamics is quit complicated, in this paper, different system designs and dynamic models will be introduced, 
simplified and more actual models and corresponding control system strategies selection and design. Ball and 
beam dynamics are shown in Figure 5. 

5.1 Modeling of basic physical sub-system model with no control involved; Ball and beam balancer 
System Modeling 

5.1.1 Ball and beam system's dynamics, basic mathematical model derivation: sliding ball. 

In this paper, different approaches for deriving the mathematical model of ball and beam will be 
introduced. Based on [1] and on system arrangement shown in Figure 2(a), dynamics are shown in Figure 5. In 
the absence of friction or other disturbances and assuming the ball is sliding, the dynamics of the ball and beam 
system can be obtained as follows; The ball moves a long the beam with acceleration, the forces that accelerate 
the ball as it rolls comes from the component of gravity that acts parallel to the beam. The equation of motion 
for the Ball is given by the following equation: 
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E F = M'ao EF = M * g * sin (a) 

Substituting and manipulating, gives: 

d"x d"x 

M * = M*g*sin(a) <=> = g * sin (a ) 

dt 2 " dr 

This model can be linearized, where for small angles, sin (a) ~~a, the model becomes: 

^— ^- = g*a (1) 
dt 

This is the basic mathematical model of the ball and beam system, this simple mathematical model shows that, 
the acceleration of the ball is proportional to the beam angle a and to gravity g. Where: x : the ball position on 
the beam, M : the ball mass, to write an approximated overall transfer function of the ball and beam system, that 
relates the input voltage and measured output acceleration, we need to replace the beam angle a, by control 
voltage, Vj m and replace the ball position x, by the position sensor output y, and finally Combine the actuator 
and sensor constants with the gravity constant, we have a single constant b , Substituting , performing Laplace 
transform and rearranging, we have the transfer function that relates the input voltage and measured output 
acceleration 

— =b *V <=> Y (s )s = b *V (s ) 

dt 2 "• " (2) 

Y (i) b Y (s ) g 



V„U) .s 2 V to (s) s 2 

Where: b: the actuator and sensor with the gravity constants. The state space representation of this approach is 
given by: 

r*,i ro nr*,i roi 

kJ"L« oJL— J*L*J" 

,r«.i 



The linear model given above is good approximation of true system dynamics, all what we need is to measure 
the constant parameter, b, it can be obtained by measured the time taken for the ball to accelerate from one end 
of the beam to the other at fixed angle of the beam, the laws of motion can be then used to integrate Eq.(l) and 
calculate constant b. The real ball and beam system has additional dynamic components due to the motor, plus 
non-linear and noise components that influence its control behavior,( coulomb friction, dead zone saturation in 
the motor input amplifier) [4]. 

5.1.2 Ball and beam system's dynamics, second approach: rolling ball without slipping. 
Based on system arrangement shown in Figure 2(a), Ball dynamics shown in Figure 5(b), and in the absence of 
friction or other disturbances and assuming the ball is rolling without slipping, the dynamics of the ball and 
beam system can be obtained as follows; 

There are two forces influencing the motion of the rolling ball, these forces are: F u : the force due to 
translational motion and F rx : the force due to ball rotation. The equation of motion for the Ball is given by: 

SF = M*a = M*g*sin(a) 

Substituting forces influencing the motion of the rolling ball , gives: 

XF = F ix +F i= M * g * sin (or ) (3) 

Now, to find forces influencing the motion of the rolling ball, and substituting in Eq.(3). Force F, x due to 
translational motion is given by: 
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F a = mx (4) 

The force due to ball rotation, F rx is found as follows: The torque produced by the ball's rotational motion is 
equal to the radius of the ball, R, multiplied by the rotational force, and is given by 

T = F rx R 

Also, the torque is equal to the ball's moment of inertia , J, multiplied by its angular acceleration, dcojdt, which 
then can be written as its moment of inertia multiplied by the double -derivative of its translational motion (x h ) 
divided by its radius, this is as shown next: 

do. 



F R 



T 

d a) b 
dt 



F R 



dt 



d (v I R) d " (x I R ) J 

J = / — = — X 

dt dt 1 R 

F ra = (J *a)/R 2 

Substituting, the moment of inertia, J, for the ball and given by ( J=(2/5)MR 2 ), we have: 

d x 



(2/5) 



M 



dt 



(5) 



Substituting Eq.(5) and Eq. (4), in Eq.(3) we have the basic mathematical model of ball and beam system that 
relates the ball displacement and angular position of the motor shaft, and given as: 



SF = F +F. = M 



: sin (a ) 



M 



(2") d x d X 

' sin (a) - — M + M 

I 5 J dt 1 dt 1 



Manipulating and rearranging, gives: 



M*g*sin(a) = dx/df[ (2/5) M + M] 
g*sin(a) = d 2 x/dt 2 [ (7/5) ] 
(5/7) *g*sin(a) = d 2 x/dt 2 

This model can be linearized, since for small values of angle a, sin(a) = a, this approximation is good for angle 
values -k/4 <a< tc/4 , sine of the angle is approximately the angle itself. Taking Laplace transform and 
rearrangingmpw? displacement, X(s), over output angular position, a(s), of the motor shaft gives : 

(5/7) ga(s)=X(s)s 2 

x CO _ (5 ii)g (6) 
a ( s ) s 

The state space representation of for this approach is given by: 

r*j_ ro r o i 

l*\\~[o oJ[^J + [ 5 /7j" 

The transfer function can be derived based on system arrangement shown in Figure 2(b), as follows; the beam 
angle controls acceleration, not position, from gravity the equation of motion for the Ball is given by: 

ZF= M*a= M*g*sin(cc) 

Assuming the ball is rolling without slipping, we have r*co= dx/dt, The moment of inertia of the solid ball is 
given by: ( J=(2/5)MR 2 ). Given centripetal acceleration with rotation of the rod as a cenlripela i = xa 2 , altogether we 



have: 



d ' x d ' a 2 d ~ x 

: x — = -g s in a 

dt dt 1 5 dt ' 



By linearization, taking Laplace transform, rearranging for transfer function, gives: 

X (s) (5 /7)g 
a (s ) s 2 



(7) 



5.1.3 The Lagrangian approach to the Ball and Beam system's dynamics: 

In the absence of friction or other disturbances, the dynamics of the ball and beam system can be obtained by 
Lagrangian method based on the energy balance of the system, The Lagrangian method is utilized to derive the 
equation of motions for the ball and beam system in the most model based research works on the ball and beam 
system [6] [11-13]. In order to derive the Euler-Lagrange equation, the first step is to define the kinetic and 
potential energy for the ball and beam. The kinetic energy of the system is given by: 



_ 
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T = T!+T 2 (8) 
Where Ti and T 2 are kinetic energies of the beam and the ball, these kinetic energies include radial and circular 
motions. The rotational kinetic energy of the beam is given by: 

T, = 0.57 ,a 2 

Where: Ji, is the moment of inertia of the beam and a is the angle velocity of the frame. The ball has kinetic 
energy, given by: 

T j = 0.5(M x 2 ) * a 2 + 0.5 Mv 2 + 0.57 ;> ft> 2 

Substituting v =Rco and 7 6 = (2/5)*M*R 2 , the rotational kinetic energy of the ball is given by: 0.5((2/5)*M*R 2 ), 
substituting, the kinetic energy of the system in Eq.( 8) will give: 

T = 0.5 ((7 + Mx 2 ) * a 2 + (7 / 5)M x 2 ) 

The potential energy of the system is exhibited by the rolling ball alone, and given by: 



M * 



■ sin or + 0.5m*g*L =i! sin£): 



P = (M*g*x+ 0.5*m * g*L)sina 

Where m : the mass of the beam ,M: the mass of the ball, L: the beam length. We can now write the 
Lagrangian, where the difference between kinetic and potential energy is called the Lagrange function, which is 
defined by (L) equation and given by: 

L=T-P 

Where : T is the kinetic energy and P is the potential energy in the system. 

L = f"o.5 ((7 + Mx 1 ) * a 2 + (7 / 5)M r 2 )1 - [(M * g * x + 0.5*m * g*L ) sin a ] 

The dynamics equation representing the variation effect of system variable is given by: 

d ( dL ~\ SL (Q) 

— |— |- =Q W 

dt { dq J dq 

(J ) d 2 x 2 (10) 

= + M I + M g * sin a - m * x * a ' v J 

U 2 ) dt 2 

Where 0: the motor output position angle, a: the beam angle, x: ball coordinate position, m: beam mass . As 
shown in Figure 2 (a) and Figure 7, It should be considered, that the beam angle, a, and DC motor output shaft 
angle, 9, are not the same, but there is also relationship between these two angles that can be seen from Figure 
6, this relation can be approximated as linear, where the distance traveled by the beam at radius equals to L, is 
equal to distance traveled by gear at radius, r, that is: 



T0 = La 



for r << L 



(ID 




Output gear angle 



DC motor & gear' 
Figure 2(a) 




Figure 7 the beam angle, a, and DC motor output position angle ,0 



Ball and beam system's dynamics in terms of ball radius R, andr : based on system arrangement shown in 
Figure 2 (a), dynamics shown in Figure 5, and in the absence of friction or other disturbances, we have: 



X — r <=> CD — + 



-+ 6 t> 



+ (x 0) 



1 ,1 2 1 •! 

U = —mv + —I a> + —I 8' 

2 2 2 



Writing and substituting the Lagrangian equation, gives: 
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if , . , \x .1 .,1 

L = — { m [x ~ + (x 9) ] + / I — + 9 I + / 6 \ 
2 I L r J * J 

d ( d \ d 

L - L = m g sin 9 

dt \ dx ) dx 

< 'J C^X 

|^ m + — — J x + ^ — — J 6 - m x = in g sin # 



Taking Laplace transform, Substituting / ( = (2 / 5)M R 2 , and manipulating for transfer function gives: 



G (s) = 



@(j) 



2 m /i " .i 
m H 



2 ( R \ S 
1 + — I — I 



The transfer function in terms of input gear angle, 0, to output ball position, X , also can be written as follows: 
Linearization of Eq.(lO) about the beam angle, a=0, gives 



= I + M I 

[r 2 dt' 



■ M g * sin a — m * x * a ~ 



( j 



\ d x 

M I = -M * g * a 

J dr 



Linearization of this equation, substituting Eq.(l 1) taking Laplace transforms, and manipulating for transfer 
function gives: 

X (s) M * g *d 1 



f J \ d X d 

+ M I = -M * g * — 9 

(R 1 ) dt 2 L 



9U) ( J \ s 

L + M 

U 1 ) 



The state space representation of this model is given below, here it is important to notice in these models, that 
by combining the ball radius R, and r constants with the gravity constant, g, we have a single constant b, this 
will give the basic transfer function model previously obtained, this is shown in below comparison of different 
approaches for model derivations: 
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The transfer function between the ball position and rotor angle can be obtained based on basic transfer function 

Y \ _ Ja_ _ A_ and on relation given by Eq.(l 1), and given by 

v ., (s ) s 1 s 2 

y ) r g r * g 1 



V in ( s) L s L s 

5.2 MATLAB/Simulink representation, simulation and analysis 

The loop system Simulink models for the simplified models is shown in Figure 8. The Simulink model 
of mathematical model given by Eq.(6) is the same as Simulink model shown in Figure 8 , but the constant b , 
now is more accurately defined and is given by (5/7)*g, the Simulink model of open loop ball and beam system 
is shown in Figure 8(c), The system response using MATLAB, can be found using the following codes 
» g=9.8; num=[5*g/7]; den=[l 0]; G_open=tf(num,den); printsys(num,den); step(num,den); 
[A,B,C,D]=tf2ss(num,den); Figure , step(A,B,C,D) 



Now , Based on [14], the nonlinear Lagrangian equation of motion is directly modeled, where Eq.(10) gives 
d(x)/dt as a function of the state and input variables , x, d(x)/dt , alpha, and d(alpha)/dt. The reference made use 
of the nonlinear function block to express this function; the function takes the input vector u=f x dx/dt a da/dt 
] and returns d 2 x/dt 2 , where each component is referred to as u[l], u[2], etc, and correspondingly, u[l]=x, 
u[2]=d(x)/dt, u[3]=alpha, and u[4]=d(alpha)/dt, the function is given by: 
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(—1/(7 I R + M ))(M g sin a 



d 'a 

M x ) 

dt 2 



And can be written in MATLAB as: 

(-l/(J/(R A 2)+m))*(m*g*sin(u[3])-m*u[l]*(u[4]) A 2) 
This model is shown in Figure 9 (a), first it is required, to define ball mass m, Ball radius, R, beam length, L, 
gear radius, r, and ball's moment of inertia / , (e.g. m= 0.1 1 1; R = 0.015; g = -9.8; L = 0.6; d = 0.03; J = 9.99e-6). 
Running all Simulink models for derived system's models, when subjected to random input, will result in the 
same position, speed and acceleration response curves shown in Figure 10(a). In Figure 10(b) is shown ball 
position and acceleration open loop response when subjected to step input. 
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Figure 8(a) Simulink model of simplified representation 
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Figure 8(c) Simulink mode of open loop ball and beam sub-system 
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Fig, 9 (a) The nonlinear Lagrangian equation of ball and beam open loop motion sub-system model 
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Beam angle VS time Ball position VS time x 10 4 Ball position VS time 
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Seconds Seconds Seconds 

(a) Response for randomly changing input (b) Response for step input 

Figure 10 open loop, ball position/time, speed/time and acceleration/time response curves 



The next MATLAB code can be used to define system's parameters in MATLAB and return transfer function 
and open loop step response 
clc, clear all 

M=input( 'Enter mass of the ball, M = '); 
m=input('Enter mass of the beam, m = '); 
R=input('Enter radius of the ball, R = '); 
d=input('Enter lever arm offset , d = '); 
L=input('Enter length of the beam , L = '); 
J=input( 'Enter half's moment of inertia , J= '); 
g= -9.8; % gravitational acceleration 

s = tf([l 0],1); G_openl = -(M*(g)*d)/(L*((J/R A 2)+ M)*s A 2); step(G_openl) 
% another simplified way of obtaining the transfer function 

K = (M*(g*d)/(L*(J/R A 2+M))); num = [-K]; den = [10 0]; G_open2 =tf(num,den); step(G_open2); [A,B,C,D]= 
tf2ss(num,den); step(A,B,C,D); 

% Obtaining closed loop transfer function with unity feedback , no gain 

G_close=feedback(G_open 1,1); axis( [ 1 , 1 00,0,2] ) ; 

title ('step response of closed loop system with H(s)=l ') 

5.3 Modeling of actuator sub-system dynamics 

The DC motor is an example of electromechanical systems with electrical and mechanical components, a 
simplified equivalent representation of DC motor's two components with gears attached are shown in Figure 
2(c) up . DC motor turns electrical energy into mechanical energy and produces the torque required to move the 
beam to the desired angular position, L , [15]. In [16] a detailed derivation of DC motor electric and mechanical 
parts dynamics are introduced. Based on the Newton's law combined with the Kirchoff s law, the mathematical 
model in the form of differential equations describing dynamics characteristics of the armature controlled 
PMDC motor can be derived, and correspondingly DC motor transfer function open loop transfer function 
without any load attached relating the input voltage, V in (s), to the motor shaft output angle, d m (s), to be given by: 

6(s) Kin 

G M „ (* ) = = 

V,„(.v) {j +R t )(J m , +b m )+ K f K s ]} (12) 

6 (s) Kin 



The PMDC motor open loop transfer function relating the input voltage, V,„(s), to the angular velocity, co(s), 
given by: 

m(.s ) Kin (13) 



,(»)=• 



v i,A s > {[(£„■/. s 2 + (RJ„ +b m L.)> + (R.b m + K : K h ) J 



Issn 2250-3005 1 1 I |September| I2013 



Page 66 



Mechatronics Design Of Ball And Beam . 



Modeling nonlinearities: Coulomb friction is a non-linear element in which forces tend to appose the motion of 
bodies in contact in mechanical systems, it acts as disturbance torque feedback for the mechanical system, 
Coulomb friction is considered to be a constant retarding force but is discontinuous over zero crossings, that is, 
when a DC motor reverses direction it must come to a stop at which point Coulomb friction drops to zero and 
then opposes the reversed direction. In effect Coulomb friction is constant when rotational velocity is not zero. 
Based on differential equations describing mechanical characteristics of DC motor; the sum of the torques must 
equal zero, we have: 

Y 1 T = J*a= J*d 2 9/dt 2 
Coulomb friction and dead zone friction, where (T ;oflrf =0) , we have: 

K t *i a = T a + T m + Tioad +Tf 
Where: i a armature current, K, : torque constant, T f : Coulomb friction torque, taking Laplace transfom: 

K*I(s) - J m *s 2 6(s)- b„*s 9(s) _T f = 
At steady state conditions, d/dt =0, gives: 

K, % = -b*co -T f 

5.4 Simplification of open loop PMDC motor system transfer functions models. 

Referring to [16] and based on the fact that, the PMDC motor response is dominated by the slow 
mechanical time constant, where the electric time constant is much faster (e.g. ten times) than the mechanical 
time constant, this can motivate us to assume that the armature inductance, L a is low compared to the armature 
resistance, R a . neglecting motor inductance, (L a =0) , will result in the following simplified first order form of 
PMDC motor transfer function in terms input voltage, V,„(s) and output speed , co,„(s) given by: 

co(s) K t 

G . ) = = 

V„0) {[(RJ m )s + (Rb m ) + K,K b ~\} 

Rearranging this first order equation into standard first order transfer function form yields: 

r < s a(s) K - K - i \\ R - J ) K " ( 14 ) 

G ( s ) = = / | .i + 1 | = 

V. (i) (Rb+K t K k ) ll [R b + if K b ) J J n + 1 



A simplified first order form of PMDC motor transfer function in terms of output angle can is also be obtained 
by substituting (L a =0), and given by: 

K , 

(15) 

= ■ = = 

V /„< J ) \\LJ s' + (RJ +bL)s'+(Rb +KK)s]\ \R J _ J 2 + (R b_ + K K )s 1 1 f KK. W s ( s + ") 

lL ' ■ ■ " " ' * " * -II L J s\s + 1 A + || 

L J , I '" R . Jj 

by substituting, (L a =0), motor equation can be simplified to first and second order system relating input voltage 
and output angle, as well as equation relating input voltage and output speed, to be given by: 

9 ) K a (s ) K , 

G , (i ) = = ' . G , (s ) = = '■ 

V ( ,(s) .v (R a J m s+K l K b ) V„(») (R.-I.n +K t K h ) 

The transfer function relating input voltage and output angle can be more simplified to have the next form: 

0(s) 1 

G „„. C* ) = = 

K,„(s) [( S J R ) 1 
j | | — - — - |+ n * K 

The geometry of the mechanical part determines the moment of inertia, the total equivalent inertia, J equ j V and 
total equivalent damping, b equ i V at the armature of the motor with gears attaches, are given by Eq.(16), The 
inertias of the gears , sensor, beam and ball have to be included in the calculations of total equivalent inertia and 
damping, for simplicity are given by Eq.(17): 

K„„ = K + b M | — | o J,,,,, = i. + i M I — I K ' 

J = J + J o / = (m/ *) /12 + (2*M * R 2 ) / 5 (17) 



5.4.1 DC motor sub-system simulation and analysis. 

The derived mathematical models of DC motor, simplified and more actual, can be represented in 
Simulink, as shown in Fig 11, with applied input voltage and output angular position, angular speed, current 
and torque; where in Figure 1 1 (a) is shown Simulink model of more accurate DC motor mathematical model 
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and the corresponding function block window is shown in 11( b). The Simulink model of simplified 
mathematical model of DC motor is shown Fig 1 1(c) and the corresponding function block window is shown in 
Figure 1 1 (d). These models can be used to analyze and evaluate DC motor sub-system performance, as well as, 
designed control system. Defining in MATLAB DC motor parameters used, and running the Simulink model 
will return angular position/time, angular speed/time, current/time and torque/time response curves shown in 
Figure 11 (e). The following code can be used to plot response curves:load Motor ,load Motor 1, load 
Motor2,load Motor3,subplot(2,2,l),plot( motor_angle ), xlabel(' Seconds'), ylabel(' Position \theta'),title(' Shaft 
angular position VS time '), grid,subplot(2,2,2),plot( motor_speed ), xlabel(' Seconds'), ylabel(' Speed \omega 
'),title(' Shaft angular Speed VS time '), grid, subplot (2,2,3),plot( motor_torque ), xlabel(' Seconds'), ylabel(' 
Torque N.m'), title(' Torque VS time '), grid, subplot(2,2,4),plot( motor_current ), xlabel(' Seconds'), ylabel(' 
Current Amp'),title(' Current VS time '), grid, 
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Figure 11 ( a) DC motor simulink model 



Figure 11 ( b) function block model 
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Figure 1 1 ( e) DC motor position/time, angular speed/time, current/time and torque/time , both for step and 

random response curves 

5.4.2 Overall Ball and beam open loop system simulation and analysis 

As shown in Figure 12(a), building overall Ball and beam open loop system Simulink model, can be 
accomplished by connecting ,in series, both open loop transfer functions of DC motor and ball and beam sub- 
systems, given by Eqs.(18) .all three proposed DC motor Simulink models, can be merged in one function block 
model shown in Figure 12(b), that can be connected to ball and beam subsystem model, and any DC motor 
model can be activated by using Simulink manual switches, based on all this, the open loop Simulink model of 
ball and beam system will have the form shown in Figure 12(c). Depending on required accuracy, selected 
control system type, algorithm and design requirement, the input signal can be connected to any of three DC 
motor models. This proposed open loop ball and beam model gives designer readings (numerical and graphical) 
about both sub-systems' dynamics, ball and beam as well as, DC motor subsystems. By defining in MATLAB, 
both DC motor and ball and system parameters, and running the Simulink model will return the following 
response curves of DC motor : angular position/time, angular speed/time, current/time and torque/time response 
curves, and the following response curves of ball system: ball position/time, ball speed/time, ball 
acceleration/time responses, applying random input signal ( to mimic ball position and corresponding shaft 
angle response) will result in all these response curves are shown in Figure 12(d) 
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Figure 12(a) Open loop simulink model of ball and beam system. 
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Figure 12 (c) Open loop Simulink model of ball and beam system 
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Figure 12 (d) angular position/time, angular speed/time, current/time and torque/time response curves of 
used DC motor and ball position/time, ball position/time, ball speed/time, ball acceleration/time 

responses, both for random input. 

5.5 Control system selection, modeling, design, and simulation 

Due to fact that it is double integrator (two poles at the origin) system and is open-loop unstable, in 
order to stabilize the system a special control system is to be designed. Ball and beam system, provides two 
outputs measurements and correspondingly two sensors one to measure ball translational position, x, and other 
to measure motor shaft angular position, 8. Correspondingly, the control system design requires two control 
loops design, inner for motor angle control, to make motor follow the input angle and other outer loop to control 
the ball position on the beam, to make ball position follow the reference input position. These controllers can be 
designed as one or separately, mainly if the motor system is required to be faster than the beam controller. 

5.5.1 Proportional Derivative, PD- controller: 

The transfer function of PD-con trailer is given by Eq.(19), Rearranging, we have the form given by 

Eq.(20): 

G PD (s) =K P +K D s (19) 

G PD (s) = K„ + K D s = K D [s + (K p IK D )] = K D (s + Z pD ) (20) 

The PD-controller is equivalent to the addition of a simple zero at Z PD =Kp/K D . The addition of zeros has the 
effect of improving the transient response 

Systems design with prefilter: Prefilter is defined as a transfer function G p (s) that filters the input signal R(s) 
prior to calculating the error signal. Adding a control system to plant, will result in the addition of poles and/or 
zeros, that will effect the response, mainly the added zero, will significantly inversely effect the response and 
should be cancelled by prefilter, therefore the required prefilter transfer function to cancel the zero is given by 
(21). 



JPefilter 



.(s) =ZJ(s+ ZJ 



(21) 



5.5.2 Inner Motor PD controller loop with deadbeat response 

Often, the goal for a control system is to achieve a fast response to a step command with minimal 
overshoot. Deadbeat response means the response that proceeds rapidly to the desired level and holds at that 
level with minimal overshoot. The characteristics of deadbeat response include; Zero steady state error, Fast 
response, (short rise time and settling time) , percent overshoot greater or equal to 0. 1 % and less or equal 2 % 
and minimal undershoot, less than ±2% error band. In controller design for deadbeat response approach, the 
controller's coefficients, (controller gains, poles and zeros), depend on the physical parameters of the system. To 
determine the optimal coefficients, that yield the optimal deadbeat response, the system's overall equivalent 
closed-loop transfer function, T(s) is compared with standard, of corresponding order, and normalized transfer 
function, (particularly the characteristic equations are compared) These coefficients were selected to achieve 
deadbeat response and minimize settling time T s and rise time, T R to 100% of the desired command. The desired 
Standard second order closed-loop transfer function for achieving desired deadbeat response specifications is 
given by Eq.(22): 
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T ( s ) = — 



(22) 



s + a co n s + co 

Based on [18], the following derivation and design for PD controller with desired deadbeat response 
specifications in terms of desired output angular position. Considering that system dynamics and disturbance 
torques depends on application shape and dimensions (robot arm, conveyer etc), the mechanical DC motor part, 
will have the form: 

K t *i a = T a + T m + Tioad +Tf 

The coulomb friction can be found at steady state, to be: ( K, *i a - b*co =7^ ). In the following calculation the 
disturbance torque, T, is all torques including coulomb friction, and given by: T=Ti oad +T f . Applying PD 
controller with deadbeat response design for output desired output angular position, the open-loop transfer 
function of the PMDC, is given by: 

K I n 

c „„„ O ) = ; 

(Lj + R„)(J„s + b m s) + (L t s + X,)(T ) + K b K t 

Manipulating for forward and closed loop transfer functions, including disturbance torques, sensor and gear 
ratio, gives: 

K nK K 
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nK K n K K K s + n K K K „ 
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n K K t K D s + n K K ( K p 
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Comparing with standard normalized third order system and manipulating, for K D and K P , given by Eq.(22) 



by: t (s 
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For design consideration, and based on required design accuracy we can use the simplified second order model, 
assuming La =0, and manipulating for closed loop with PD controller, gives: 
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G (s ) = 

V ,„ (s ) J [(if )i + (J? „&,,„,„ + X , K „ +iiJf K f K D )] + nK fn K : K p 

Comparing with standard normalized second order system and manipulating, for K D and K P , gives 

nK K „ + nK K ., s 



r i r KK b \\ 

s | j + 1 b + | 

L I " R ., )\ 

K t K r + K t K D s 

V Ui (s) s [(R t J M )s + (R J) m +K t K b )] 



(RJ ) 



^ R «*«,»(, +K , K I, + " K ,,„, K , K D ^ nK „ K , K 1 

S + I — I s + 



5 + aco s + co 



R J 

a eqwv 

R,b + K K b + nK iui KK i 
R J .. . 



( R J«„„ r ) 
ae) n RJ tqu . y -(R,b +K : K b ) 
nK K 
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nK K K „ R J 0)' 

2 pot iP a tqim n 

co _ = o K = 



(R J . ) nK K 

v a equiv ' pot I 

z 

PD 

Referring to [18] The controller gains K p and K D depend on the physical parameters of the actuator drives, to 
determine gains that yield optimal deadbeat response, the overall closed loop transfer function T(s) is compared 
with standard second order transfer function given by (22), and knowing that parameters a and co n are known 
coefficients of system with deadbeat response given by [13], a = 1.82 and co n T n = 4.82 , T n =2 and gives the 
following: co n T n = 4.82 , co n = 4.82/2=2.41. In [18] to simplify and accelerate Mechatronics design process, a new 
MATLAB built-in function with specific purpose, which is design and verification for desired deadbeat 
response specification, particularly for desired settling time. The presented built-in function is named 
deadbeat(a), the input argument a, can have the value of 1 or 2. The input argument 2 for controlling electric 
motor system output angular displacement with PD-controller design for deadbeat response with desired settling 
time, this new built in function determines the coefficients (gains pole and zero) that yield the optimal deadbeat 
response for desired output speed or angle control, as well as performance specifications in terms of Mp, E ss , Ts, 
final output value. 

5.5.3 Inner Motor PD Position control loop with position and rate feedback 

Considering the output position control, applying position and rate feedback ,( see Figure 13(b)) , to 
have the motor follow the desired output angle, the input voltage to DC motor will be given by Eq.(23). 
Rearranging, we have Eq.(24): 

V, = K p [9 t -9 ]-K D a> » V, = K p [9, - 9 ] - K D s 9 c (23) 

(24) 

Finding the overall equivalent closed loop transfer function, and comparing with corresponding standard 
transfer function ( second or third order), and by comparison, we can find the gains K P and K D depending on 
motor parameters and desired damping ratio, C, and settling time, T s . Standard second order system, in terms of £ 
and <d„, is given by Eq.(25): 

lh_ = a " (25) 



T (i ) = = — " - 

9 ialred s " + 2£<u n + o? ~ 

Assuming the constrains for this system are to be , desired damping ratio of 0.707, or overshoot and Peak time 
of T P =0. 1 second, we find undamped natural frequency co„, and correspondingly K P , K D 

- ln(%OS / 100) 7i 
■^n 1 + \n 2 (%OS /100) m n ^\ - f 2 

From, the DC motor transfer function, depending on desired accuracy, and correspondingly applied 
mathematical model, and by substituting the value of V-„ we have: 

\(sJ R ) 1 
V . = s | | - - \ + n*K t | O ) 

IA B ***J J 

\( S J R ~) 1 

V. = K p - (K r - K D s )0 o 0,(i)s| \ + n * K t \ = K p (s ) - ( K p + K D s)6 (s ) 

1{"* K >J J 

\( S J R \ 1 

e (s )j 1 1 - " \ + n * k : \ + (k r + k n s)g i u) = k r e i u) 

IK H * K > J J 

( \~fsJR \ 1 1 

6 (» ) s | | \ + n*K, | + {K , + K s ) = K r 9. (s ) 



OAs) ( s 2 j R } 
" — -+ (nK , + K )s + K r 

By comparison with standard second order system, we find PD controller gains K P , K D 



5.5.4 Inner Motor PD Position control with only position feedback 
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Considering the output position control, applying only position feedback with proportional and 
derivative terms, as shown in Fig 13(a). The input voltage to DC motor will given by 

To find the gains K P and K D , apply the same procedure, for desired damping ratio, C, and settling time, T s . 

5.5.5 Inner Motor PID Position control loop with only position feedback 

Considering the output position control, applying only position feedback with proportional, derivative 
and integral terms, as shown in Fig 13(b). The input voltage to motor will given by: 

k \e - 9 l 

V, = K r [ff. -g a ] + S K D + ■ — 

S 

The gains K P , K D and K b are found by applying the same procedure, for desired damping ratio, £ and settling 
time, T s . 
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Figure 13(a) Motor PD Position control with position and/or rate feedback ( to switch use manual switch) 
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Figure 1 3 (b) Motor PID Position control with position feedback 
5.5.6 Outer Ball and beam PD, PID Position control with only position feedback 

The same procedure can be applied, to calculate controller gains, for controlling beam angle, when applying PD 
or PID, this is shown in Fig 13(c) 



5.5.6.1 Outer ball and beam PD Position control loop with position and rate feedback 
The input beam angle, a, to ball and beam system, applying PD controller, will be given by: 



a = K 



dx 



dt 



Taking Laplace transform, 

a = K p [X ,(»)-* ,(»)]-«,<(!) 

Substituting a, ball and beam transfer function, as shown below, we have 



X (>) (5/7)f 



a (s ) 



•» a ( s ) 



X (s) <=>«(*) = — X (») 

(5)(9.8) 7 



-X „ (.v ) = K „ [X . (j ) - X „ (s )]-f„j(s) 



IK 



X . (i) i +7*f a s + 7*-, 

This is the overall equivalent closed loop transfer function of ball and beam, comparing it with standard second 
order transfer function, and by comparison, we can find the gains K P and K D depending on plant parameters and 
desired damping ratio, t, and settling time, T s . as shown next: 



t co 



X (O 



X (s ) s + 2^co + co 



IK 



+ IK „ j +1K , 
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K 



7 7 

Assuming the constrains for this system are to be , desired damping ratio of 0.707, or overshoot and Peak time 
of T P =1 second, we find undamped natural frequency co„, and correspondingly K P , K D 

- ln(%OS / 100) 



Solving we have: co„ = 4.4422, k 



^tr 1 + \n 2 (%OS /100) 
2.8190 ^ K = 0.89728 




Figure 13 (c) Ball and beam position control with position feedback 



5.5.7 Sensor sub-system modeling 

There is a relationship between the beam angle, a, and DC motor output position angle , 8, where the 
distance traveled by the beam at radius equals to L, is equal to distance traveled by gear at radius, r, that is: 
!0 = La o a = (r i L)6 , Where: The beam angle a = (r / l)0 , and the motor angle e = (L / r )a -The 
transfer function between the ball position and rotor angle can be obtained based on basic transfer function 
given by Eq.(2) and on relation given by (1 1), and given as shown next: 

y ( s ) b s y (s ) f g r * & i 

V h} (s ) s 2 s 1 L s 2 L s 2 

The position of the ball on the beam can be measured using special sensors, e.g. linear Potentiometer, in [19] a 
linear resistive wire is used as a position sensor and the conductive ball moving along it acting as a voltage 
divider Linear resistive wire has a relative large resistance of about 2.15 ohm/m. A series resistor is added to the 
circuit in order to divide the voltage to avoid unwanted heat. According to the position (X) of the ball, the 
voltage will linearly vary from 0.1 V to 0.74 V. Schematic of the resistive wire circuit is shown in Figure 14. 




Figure 14 Schematic of the linear resistive wire circuit "position sensor" [19] 

(Voltage change) (0.74 - 0.1) (0.64) 



(Position change) 

For voltage range of 10 volts; k _ =10/40= 0.25 



( L ) 



( 40 ) 



5.6 Electronics and interfaces simulation and testing: simulation in Proteus; To test and evaluate the 
selection and integration of circuit design, programming and interface components, simulation using ISIS- 
Professional Proteus is used. The control program written in C language, with the help of MikroC program is 
converted to Hex. File and downloaded on the simulated PIC -microcontroller and circuit, the simulation is 
shown in Figure 18, after testing, evaluating and optimizing various aspects, the final simulation results show 
the correctness of written program, interfaces, and microcontroller, all these can be used to build the physical 
circuit. 



5.7 Ball and beam overall system simulation and analysis. 

Based on the specification of requirements and design, the subsystems models and the whole system 
model, are to be tested and analyzed,. If the specification are not satisfied,, modifications can be made, if the 
specifications are satisfied the model can be Optimized. 

The primary basic Simulink model, with PID control algorithm for both loops is shown in Figure 15, 
these two controllers will be replaced with corresponding proposed control algorithm to achieve optimal 



Mechatronics Design Of Ball And Beam . 



response. Switching model, using manual switch, to nonlinear PMDC motor model, applying two Simulink 
built-in PD controllers for inner and outer loops, subjecting system to step input, and applying Simulink block 
tuning; will result in ball position/time, speed/time, ball acceleration/time ,beam angle/time responses curves 
shown in Figure 16(a)) , when system is subjected to random input, will result in ball and beam response curves 
shown in Figure 16(b). Switching model to simplified PMDC motor model, and subjecting the system to step 
input, will result in response curves shown in Figure 16(c)), comparing response curves, show that the 
overshoot using simplified motor model is less , as well as setting time, but generally system takes 10 second 
response to reach desired position which is large for such system.Applying Outer Ball and beam PD Position 
control with only position feedback, Figure 17(a) Switching model to simplified PMDC motor model, defining 
parameters, running the model, will result in response ball position/time, ball speed/time, ball acceleration/time 
and beam angle/time response curves response curves shown in 17(b). Refining the PD Position control design, 
for both loops, will result in more acceptable response curves response curves shown in 17(c), the system 
reaches desired output at 2.3 seconds. 
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Figure 15 The primary Simulink model, with PID controller. 
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Figure 17Ball and beam system step and random input response, using actual PMDC model. 
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Figure 16 (c)) Ball and beam system step response, using simplified PMDC model 
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Figure 16 (a) Ball and beam PD Position control with only position feedback 
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Figure 17(b) Ball and beam system step response 
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VI. PROTOTYPING, TESTING, EVALUATION AND OPTIMIZATION. 

There is no single model which can ever flawlessly reproduce reality, there will always be errors called 
as unmodeled errors between behavior of a product model and the actual product. In order to take into account 
the unmodeled errors and enhance precision, performance in the design process, the Mechatronics design 
approach includes prototyping phase. Prototyping development may be carried out in the following two forms; 
Virtual Prototype and Physical Prototype, in this paper, only virtual prototype is built, Virtual Prototype is 
computer model of a product presented in a virtual environment with, ideally, all information and properties 
included, to visualize text, analyze and evaluate machine movement and logical operations. Different sources 
introduce different virtual prototypes for evaluating the overall system selection, integration and design, in [22], 
introduced model shown in Figure 18(a), the software used is MATLAB/Simulink. 




Figure 18(a) overall system virtual prototype^ Figure 18 Electronics simulation in Proteus software 

VII. CONCLUSION & FUTURE WORK 

This paper proposes the conception and development of ball and beam system based on Mechatronics 
design approach, where realization of design is achieved applying simultaneous, parallel, design approach and 
consideration. A complete system's components' selection, design and integration, as well as, modeling, 
simulation and analysis are presented. The simulation results for PD controller design are shown. Depending on 
desired accuracy, different mathematical and Simulink model forms of the DC motor and plant have been 
developed. The proposed Mechatronics design and models are intended for research purposes, as well as, for 
application in educational process. The proposed models were created and verified using MATLAB Simulink, 
Proteus software. 

This paper is part (I) of authors' research study on Mechatronics design of ball and beam system, as future work 
of this research, author intends to built and test physical prototype of this design, and to propose a general 
MATLAB/Simulink model with its function block parameters window, that can be used to control ball and 
beam system using different control strategies, as well as to simplify and accelerate the selection, design and 
verification processes 
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